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I identify the template organizing the chaotic dynamics 
of a bouncing ball system. I further show how to estimate 
the topological parameter values of the system directly from 
a time series — a process I call "topological time series analy- 
sis." Two distinct methods to determine the topological pa- 
rameters are illustrated and compared — the "pruning front" 
procedure and a "braid analysis." Both procedures lead to 
compatible results. 

05.45.+b, 47.52.+J 



I. INTRODUCTION 

Braids arise as periodic orbits in dynamical systems 
modeled by three-dimensional flows [1-4]. The existence 
of a single periodic orbit of a dynamical system can imply 
the coexistence of many other periodic orbits [5-8]. The 
most well-known example of this phenomenon occurs in 
the field of one-dimensional dynamics and is described 
by Sarkovskii's Theorem [9]. Less well-known is the fact 
that analogous results hold for two-dimensional systems 
[3]. In one-dimensional dynamics it is useful to study 
the period (or the permutation) of an orbit [10]. In two- 
dimensional systems it is useful to study the braid type of 
an orbit [2]. Given this specification, we can ask whether 
or not the existence of a given braid (periodic orbit) forces 
the existence of another; as in the one-dimensional case, 
algorithms have recently been developed for answering 
this question [11-13]. 

As originally observed by Auerbach et. al., unstable 
periodic orbits are available in abundance from a single 
chaotic time series using the method of close recurrence 
[8,14,15]. By a "braid analysis" I propose to analyze 
a chaotic time series by first extracting an (incomplete) 
spectrum of periodic orbits, and second ordering the ex- 
tracted orbits according to their orbit forcing relation- 
ship. As shown in this paper, it is often possible to find a 
single periodic orbit, or a small collection of orbits, which 
forces many orbits in the observed spectrum. These or- 
bits also force additional orbits of arbitrarily high pe- 
riod. This analysis is restricted to "low-dimensional" 
flows (roughly, flows which can be modeled by systems 
with one unstable Lyapunov exponent), however it has a 
strong predictive capability. 

I would also like to point out that this analysis gives 
us an effective and mathematically well defined "pruning 
procedure" for chaotic two-dimensional diffeomorphisms 



[16]. Instead of asking for rules describing which orbits 
are missing (pruned), I instead look for those orbits which 
must be present. For low period orbits (say up to period 
11) this procedure can predict all those orbits which must 
be present in the flow. This procedure will usually miss 
orbits of higher period, however from an experimental 
viewpoint the low period orbits are the most important 
and accessible. Orbits of low period often force an infinity 
of other orbits. This is illustrated in one-dimensional 
dynamics by the famous statement "period three implies 
chaos" [17]. An analogous statement in two-dimensional 
dynamics is that a non-well-ordered period three braid 
implies chaos [18]. 

This paper is organized as follows. Section II reviews 
the dynamics of the bouncing ball system. In section III 
I identify the template organizing the chaotic flow of the 
bouncing ball system. It is a horseshoe with a full twist. 
In section IV I show how braid analysis works by apply- 
ing it to times series data generated from the bouncing 
ball system. The analysis builds directly on the origi- 
nal topological analysis of such data sets due to Mindlin 
et. al. [19]. This section also illustrates how easily mea- 
sured braid invariants of the periodic orbits lead to strong 
dynamical information about the flow — without the need 
for a problematic and detailed symbolic description of 
the orbits in phase space. Section V illustrates how a 
pruning front can be estimated from a collection of low 
period orbits. The results are compatible with the braid 
analysis of the previous section. Section VI offers some 
concluding remarks. 

In the example studied in this paper I do have good 
control of the symbolics. In principle, though, a braid 
analysis does not require good control of the symbolics 
(a good partition) and can thus overcome some of the 
current difficulties associated with finding good symbolic 
descriptions for (nonhyperbolic) strange attractors [20]. 



II. BOUNCING BALL SYSTEM 

Consider the motion of a ball bouncing on a period- 
ically vibrating table. This bouncing ball system arises 
quite naturally as a model problem in several engineer- 
ing applications: examples include the generation and 
control of noise in machinery such as jackhammers, the 
transportation and separation of granular solids such 
as rice, and the transportation of components in au- 
tomatic assembly devices which commonly employ os- 
cillating tracks. Several researchers have studied one- 



dimensional models of the bouncing ball system which 
include the coefficient of restitution (0 < a < 1), and 
many have also noted the existence of a large class pe- 
riodic, chaotic, and eventually periodic orbits known as 
"sticking solutions" [21]. More details can be found in 
Ref. [22]. All these models have been termed the "exact" 
one-dimensional model of the bouncing ball system. The 
phrase "one-dimensional" refers to the number of degrees 
of freedom the ball moves in and not to the dimension of 
the phase space model. 

To fix a notation which allows an easier comparison 
with experiments, recall that the dynamics of the bounc- 
ing ball system can be found by solving the (implicit) 
nonlinear coupled algebraic equations known as the phase 
map, 
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where 6 k = uit -\- Oq, and v k are the phase and velocity of 
the k-th impact between the ball and oscillating table, A 
and u) are the table's amplitude and angular frequency, 
a is the coefficient of restitution, and g is the gravita- 
tional acceleration. The table's forcing period is denoted 
by T = 2-k/uj. The implicit phase map and explicit ve- 
locity map constitute the exact model of the bouncing 
ball system. Earlier experimental studies showed an ex- 
cellent correspondence between the exact model and the 
dynamics of an experimental bouncing ball system, all 
the major bifurcations predicted by the exact model oc- 
curred within the experimental system [23]. Observations 
between the model and experiment agreed to within %2 
with no fitted parameters. A public domain program, the 
Bouncing Ball Simulation System, is available to simu- 
late the exact model [22]. I use this program to obtain 
the time series data analyzed here. 

Experiments illustrating chaos in the bouncing ball 
system usually proceed along the following lines. The 
amplitude of the table driving the ball is slowly in- 
creased while monitoring the dynamics of the bouncing 
ball through an experimental impact map, which is sim- 
ilar to a next return map. In essence, an experimental 
bifurcation diagram is created. The coefficient of resti- 
tution can be changed from around 0.2 to 0.8 by using 
different materials for the ball (eg., wood, plastic, steel). 
Experimentally, it is observed that a chaotic invariant set 
is seen at the end of the period doubling cascade, but for 
a further increase in the driving amplitude, the strange 
attractor is destroyed by a crisis [24]. The dynamics 



of the ball after this crisis can result in motion which 
can quickly approach a periodic sticking solution (gen- 
erally speaking, for smaller values of a), or can exhibit 
long transients — sometimes called 'transient chaos' [25] — 
following the "shadow of the strange attractor" (gener- 
ally speaking, for larger values of a). 

Direct simulation of the exact model exhibits a similar 
behavior. Figure 1 presents a bifurcation diagram show- 
ing a period doubling route to chaos for a = 0.5. Note 
that this strange attractor is approached in the same way 
as it would be in an experiment, namely, by slowly scan- 
ning the amplitude until the end of the period doubling 
cascade is reached and a non-periodic orbit is observed. 
In simulations (A = 0.012) the strange attractor is found 
to be stable for over 10 6 impacts. For A > 0.0.118 
and A < 0.019 a first crisis occurs which expands the 
size of the strange attractor. Between A = 0.0121 and 
A = 0.0122 a second crisis occurs which destroys this 
strange attractor. For A > 0.0122 the orbit follows the 
shadow of the strange attractor for a number of impacts 
but eventually converges to a sticking solution (typically 
after 10 2 to 10 3 impacts). In both experiments and sim- 
ulations, the pre-crisis (chaotic) dynamics and post-crisis 
(eventually periodic) dynamics are usually easy to distin- 
guish because the range of impact phases explored by the 
ball suddenly widens after the crisis. In the simulation 
shown in Figure 1, the chaotic dynamics is confined to a 
phase between —0.1 < 6/2-k < 0.3 where as the (second) 
post-crisis dynamics explores almost the entire range of 
phases available. This feature provides a nice signature 
to distinguish the pre- and post-crisis dynamics in both 
experiments and simulations. 

This general scenario of period-doubling, chaos, crisis, 
and sticking solutions (possibly with transient chaos) is 
not confined to a few selected parameter values but is 
generally observed for a wide range of a. Both exam- 
ples of crisis occurring in the bouncing ball system, and 
the existence of multiple coexisting attractors, can not 
be explained by a one dimensional unimodal theory, and 
provide the first indication of the need for a two dimen- 
sional theory. 



III. BOUNCING BALL TEMPLATE 

The first step in analyzing the structure of the chaotic 
set in the bouncing ball system is the identification of 
a template which captures the organization of its pe- 
riodic orbits [19]. The template is a nice (canonical) 
representation of the stretching, folding, and twisting 
of phase space resulting in a particular chaotic form. 
To visualize the template arising in the bouncing ball 
system I plot a chaotic orbit in the three dimensional 
space, (sin(o;i), v(t), x(t)), where the first coordinate is 
the table's (normalized) time dependent forcing ampli- 
tude, and the remaining coordinates are the ball's veloc- 
ity and height. 



Inspection of Fig. 2(a) reveals a band like structure 
with a half twist occurring where the ball reverses ve- 
locity when it hits the table, and an additional smaller 
fold on the outer edge of the band near the top of the 
figure. A schematic of the sheet like structure is pre- 
sented in Fig. 2(b). A template is nothing more or less 
than this sheeted structure collapsed to a single sheet 
and moved by a sequence of isotopies to a standard form. 
This sheeted structure is perhaps easier to see in Fig. 3. 
Here the pre-image of the fold can be traced back to its 
impact point with the table. The impact phase of the 
fold point is in the vicinity of 6 k ~ 0. The folding of 
the strange attractor in phase space occurs in the region 
where the table's impact velocity is maximal: roughly, 
orbits hitting at phases greater or lesser than this value 
get less of a kick from the table and hence do not travel 
as high. 

Figure 4 shows how this sheeted structure can be put 
into a template of standard form. In Fig. 4(a) the evo- 
lution of a small section of an unstable manifold (rep- 
resented by an arrow) is shown as it is carried by the 
template. To reach a canonical form I first pull the fold 
point all the way to the left (thus going from a pruned to 
unpruned system), and second identify and cut through 
the trajectory of the fold point. In the language of Cvi- 
tanovic et. al. [5], this fold point is a primary tangency. 
As shown in Fig. 4(b) each branch of the template is now 
given a symbolic label. I also put the insertion layer of 
the template in standard form (back to front) and slide 
all the branch twists to the top of the diagram [4]. The 
template of the bouncing ball system in standard form is 
shown in Fig. 4(c). At this point I notice that by sub- 
tracting a full twist from the entire template I arrive at a 
horseshoe template in standard form (Fig. 4(d)) — thus, 
the template in the bouncing ball system for the parame- 
ter range considered is the horseshoe with a global torsion 
of —1. In the next section I verify that the template iden- 
tified is correct by comparing topological invariants cal- 
culated from the horseshoe template and those extracted 
directly from a chaotic time series. 



IV. BRAID ANALYSIS 

A braid analysis of a low dimensional chaotic time 
series consists of four steps once an appropriate three- 
dimensional space is created [19]: (i) the periodic orbits 
are extracted by the method of close recurrence [22,26], 
(ii) the braid type of each periodic orbit is identified and 
the orbits are ordered by their two-dimensional forcing 
relationship [11,27] (iii) a subset of braids are selected 
which have maximal forcing and which force the orbits 
extracted in step (i), and (iv) if possible, an attempt is 
made to verify that some of the predicted orbits (not 
originally extracted in step (i)) are found in the system. 

In practice, steps (i) and (ii) are greatly simplified if 
the template or knot-holder organizing the flow can be 



identified using the procedure described by Mindlin and 
co-workers [4,19,28,29]. Knowledge of the template helps 
in obtaining the symbolic names of the periodic orbits 
and in calculating the forcing relationship for the specific 
braids in that template. For instance, if the template is 
identified as a two-branch horseshoe knot holder (as is 
the example studied in this paper), then the theory of 
quasi-one-dimensional (qod) orbits of Hall [6,27] can be 
applied to simplify the analysis. 

Although template identification is very valuable, it is 
not essential for a braid analysis. Nor is the symbolic 
identification of the extracted orbits. In the worst case 
a braid analysis does require that the the braid conju- 
gacy class of each extracted periodic orbit is identified 
(see Elrifai and Morton [30], or Jaquemard [31] for algo- 
rithms), and that the minimal Markov model (a 'train 
track' in the language of Thurston) can be constructed 
for each braid (see Bestvina and Handel [11], Los [12], 
and Franks and Misiurewicz [13] for algorithms). Algo- 
rithms exist for both of these steps, although the most 
computationally efficient version of the braid conjugacy 
algorithm is probably not an effective solution for braids 
beyond period 8. 

To illustrate braid analysis I simulated the bounc- 
ing ball system for 10 5 impacts with system parameters 
a = 0.5, u) = 2tt60, and A = 0.01215. The resulting next 
impact map, (cj) k , (j) k+1 ), (j) k = 9 k /2ir, is shown in Fig. 
5. The many leaves of this return map once again indi- 
cates that the symbolic dynamics of this system should 
exhibit departures from that predicted by a one dimen- 
sional unimodal theory. The inset of the Fig. 5 shows an 
expand view of the region surrounding the maximum of 
the map. Three distinct leaves are visible in this region 
and this suggests that, to a first approximation, the sym- 
bolic dynamics of the system should be describable by a 
a three step pruning front. 

To extract the (approximate) periodic orbits by the 
method of close recurrence I first convert the next im- 
pact map from the sequence of values (<f> k , <f> k +i) directly 
into a symbol sequence of 0's and l's. In this particular 
instance, I found that an adequate symbolic description 
(at least up to period 11 orbits, or approximately one 
part in 2 11 ) is obtained by choosing the maximum value 
of the next impact map at the three leaves shown in Fig. 
5. Orbits passing to the left of the maximum in the vicin- 
ity of a given layer are labeled zero, and those to the right 
are labeled one. Next I search this symbolic encoding for 
each and every periodic symbol string. Every time a pe- 
riodic symbol string is found I calculate its (normalized) 
recurrence and then save the instance of the orbit with 
the best recurrence. The advantage of this procedure of 
orbit extraction is that it is exhaustive. I search for ev- 
ery possible orbit up to a given period. In these studies 
I searched for all orbits between periods 1 through 11. 
Some of the extracted periodic orbits are shown in Fig. 
6. 

The resulting spectrum of periodic orbits up to pe- 
riod 11 is shown in Table I. Simple topological invariants 



(linking numbers and relative rotation rates) of the ex- 
tracted orbits are calculated and compared with those of 
a horseshoe with a global torsion of -1. There are no dis- 
crepancies. This indicates that — at least to this level of 
resolution — the template is correctly identified and the 
symbolic partition is adequate. The orbits which are 
present in (the full shift) complete hyperbolic system, 
and not present in the Tables in Appendix A, are said to 
be ■pruned. Our goal is to predict as best as possible the 
pruned spectrum from the chaotic time series. 

The symbolic label (up to braid type) can also be de- 
termined by considering simple and easily computable 
braid invariants. For instance, as pointed out by Hall 
[27], the exponent sum (simply the sum of braid cross- 
ings in the example) is a complete invariant for horseshoe 
braids up to period 8 (see Table II). Also, as inspection 
of the tables in Appendix A reveals, the exponent sum 
manages to distinguish most of the the pseudo-Anosov 
horseshoe braids of periods 9, 10, and 11 as well. Thus, 
I see that an easily determined quantity measured from 
a time series leads to the conclusion that the flow con- 
tains a chaotic invariant set — with out the calculation 
of more detailed quantities such as fractal dimensions or 
Lyapunov exponents. 

The goal of a braid analysis is to find a small subset 
of orbits, called a "basis set" [7], which forces the the 
observed periodic orbit spectrum. One sensible way to 
proceed in identifying such a collection of orbits is to 
calculate the spectrum of orbits forced by a few high 
entropy orbits to see if they can capture most or all of 
the observed spectrum. If some orbits are left out then 
they are systematically added to the basis set until all 
the orbits in the observed spectrum are captured. 

Using the tables in Appendix A I find that the high- 
est entropy orbit in this particular data set is Sg (/12 = 
0.397) which happens to be a quasi-one-dimensional or- 
bit. Thus, using the results of Hall [6,27], I calculate 
the forced spectrum of this orbit by one-dimensional uni- 
modal kneading theory [32]. I find that up to period 
eleven the s| braids forces s^ x (s^), s\, s^ x , Sg, s\ x , s^ , 
s§, s^Q, Sq, s\ , Sg, s\, s\, s\. There are still many or- 
bits unaccounted for in the observed spectrum. I thus 
examine the orbit with the next highest entropy. It is 
the s| (h 2 = 0.377) orbit. The s| braid is not QOD. Its 
forced spectrum can be calculated either by obtaining 
the train track for this braid (and the associated Markov 
model) [11], or by the method of pruning fronts briefly 
illustrated in section V. Using a train track calculation 
I find that the s| braid forces s^ (s^), sf (s\ ), Sg, 

s ii ( s ii)> 4, 4, *ii. s io> 4, 4, 4, 4- Comparing the 
spectrum forced by the union of these two orbits with 
the observed spectrum (Table I), I find that only one or- 
bit is unaccounted for, the finite order braid s|, which is 
the maximal orbit in the observed data set in terms of 
one-dimensional unimodal theory. 

Adding this orbit to the collection I determine that a 
basis set which accounts for the observed spectrum — up 
to braid type — is {s|, s|, s|} 



V. PRUNING FRONTS FROM PERIODIC 
ORBITS 

I now attempt to predict the (low period) forced orbits 
by using the extracted periodic orbits to systematically 
construct an approximation to the pruning front [5]. The 
braid analysis of the previous section only specifies the 
existence of orbits up to braid type. Thus, for instance, 
it might only predict the existence of one individual pe- 
riodic orbit in a given saddle-node pair. The pruning 
front procedure is more specific, it actually forces indi- 
vidual periodic orbits as denoted by their complete sym- 
bolic label. Not unexpectedly, therefore, the "basis set" 
of periodic orbits needed to constructed an approximate 
pruning front may be larger than that found in a braid 
analysis. 

As a first step in obtaining an approximate pruning 
front I plot the trajectory of a single chaotic orbit in the 
symbol plane [5]. The data is a symbolic symbol string 
constructed in the previous section of the form 

S = . . .S_3S_2S_lS .SlS2S 3 . . . 

where symbols to the left and right of so are the past 
and future symbols respectively. The coordinates of the 
symbol plane for a horseshoe are calculated from the well 
ordered past (c^) and future (bi) symbols as follows: 
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If s is an infinite symbol string generated by a chaotic 
orbit, then D is infinity in the above sums. However, 
since I am dealing with finite data sets, I approximate the 
symbol plane coordinates of a point by taking D = 16. 
In this way I can use a finite symbol string from a chaotic 
trajectory to generate a sequence of points on the symbol 
plane. The resulting plot for the data is shown in Fig. 7. 

An expanded view of the primary pruning front region 
(center right of full diagram) is shown in the inset of Fig 
7. The two dimensional nature of the data set is indicated 
the by the steps in the diagram. If the data set was one 
dimensional then a vertical pruning front with no steps 
would be seen. Such vertical pruning fronts are found, for 
instance, when the bouncing ball system is much more 
dissipative (a ss 0.1). As a rule of thumb, the depth of 
the steps increases as the dissipation decreases. In this 
example, the steps are easier to see in the iterates of the 
pruning front. 

Now, to construct an approximate pruning front I plot 
all the periodic orbits (periods one through eleven) ex- 
tracted in the previous section on the symbol plane and 



examine their location in the region of the pruning front. 
This is shown in Fig. 8 in the same region as that found 
in the inset of Fig. 7. The periodic orbits closest to the 
right and the center (i.e. closest to the pruning front 
suggested by Fig. 7) are selected as a basis set for con- 
structing an approximate pruning front. Labels for these 
innermost periodic orbits are indicated in Fig. 8, and the 
last digit in the symbolic label (for the saddle-node part- 
ners) is determined by whether rightmost point of the 
orbit lands above (1) or below (0) a line through the cen- 
ter of the symbol plane. To construct the approximate 
pruning front I take the orbit in each saddle-node pair 
which is larger (rightmost) by unimodal ordering. Thus 
in this example I construct the pruning front from the set 
of periodic orbits (from bottom right to center): s|(0), 

4,(0), ^(o), s |(o), 4(1). 

An approximate pruning front is then constructed from 
a continuous sequence of horizontal and vertical line seg- 
ments connecting these periodic orbits: 
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(011101101.101101110,011101101.101111011), 

V4[s| ,sj] -> 

(011101101.101111011,110111101.101111011). 

By construction, this pruning front generates the same 
periodic orbit spectrum — up to period eleven — as that 
recorded in Table 1. 

Like the braid analysis, beyond period eleven this 
pruning front should begin to generated fewer orbits than 
are actually present. Both of the pruning methods illus- 
trated here are systematic approximations in the sense 
that, given a periodic orbit spectrum up to period P, 
these methods generate an exact spectrum up to some 
period Q, beyond which both methods then provide lower 
bounds on the periodic orbit spectrum. 



VI. CONCLUSION 



a dissipative bouncing ball system. Two distinct tech- 
niques where used to predict orbit forcing — a braid anal- 
ysis, and the pruning front approach. Both techniques 
provide an effective procedure for calculating the orbit 
spectrum of low period orbits. Both procedures also pro- 
vide information (a lower bound) for the spectrum of all 
periodic orbits. 

Each procedure for estimating the periodic spectrum 
has distinct advantages and disadvantages. The braid 
analysis does not require a symbolic partition, and is 
thus useful in the cases (eg. low dissipation) where de- 
termining an exact, or approximate symbolic partition, 
is problematic. The braid analysis is also based on a 
rigorous mathematical foundation. The braid analysis, 
though, only provides information about orbit forcing up 
to braid type. The chief advantage of the pruning front 
approach is that it provides information about individual 
periodic orbits. Its chief disadvantage is that it, as yet, 
rests on a weaker mathematical foundation and requires 
the construction of a symbolic partition. 

In retrospect, I find it remarkable that such a small 
subset of periodic orbits (which are rather easy to get 
from experiments) contain so much topological and dy- 
namical information about a (low-dimensional) flow. A 
few low period orbits are sufficient to determine the tem- 
plate describing the stretching and folding of the strange 
set. The template provides an upper bound to the topo- 
logical entropy and is, in a sense, a maximally (i.e., a full 
shift) hyperbolic set which can be formally associated to 
a (possibly nonhyperbolic) strange set. In this paper I 
show how periodic orbits (and their associated hyperbolic 
sets) can be used to obtain an approximation to a strange 
set which is probably not hyperbolic. Formally, I might 
say that the hyperbolic set associated with each pseudo- 
Anosov braid is embedded within the strange attractor 
I am trying to describe in the sense that the (possibly 
nonhyperbolic) strange set must contain at least all the 
orbits forced by the extracted pseudo-Anosov braid. 
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I have illustrated how to determine the topological 
form (the template) and how to estimate topological pa- 
rameters directly from a chaotic time series generated by 



APPENDIX A: INVARIANTS 

Topological invariants of horseshoe braids from period 
one through eleven are shown. 
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FIG. 1. Bifurcation diagram for the bouncing ball system. 



FIG. 2. (a) Three dimensional plot of a chaotic trajectory, 
(b) Schematic of the sheeted structure governing the evolution 
of the chaotic trajectory. 



FIG. 3. Two dimensional plot of a chaotic trajectory in 
the bouncing ball system. Inset shows template governing 
the evolution of the orbits. 



FIG. 4. Transformations taking the template found in the 
phase space to a template in standard form. 



FIG. 5. Next phase map for the bouncing ball system. In- 
set shows expanded view of the region near the maximum. 



FIG. 6. Periodic orbits extracted from a chaotic time series. 
The exponent sum identifies the orbit up to braid type. The 
linking number of the Z(0, 10) orbits is also indicated. 



FIG. 7. Symbol plane diagram generated by a chaotic tra- 
jectory in the bouncing ball system. Inset suggests an outline 
for a pruning front. 



FIG. 8. Extracted periodic orbits, up to period eleven, 
plotted on the symbol plane. Taking the innermost points 
leads to an approximation of the pruning front. 



TABLE I. Spectrum of low period orbits extracted from a 
chaotic time series of the bouncing ball system (all orbits with 
e < 0.01) are shown). Extracted orbits and their (best) nor- 
malized recurrence are recorded. Note that in this particular 
example all saddle-node partners are detected in pairs. 
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TABLE II. Exponent sums for horseshoe braids up to pe- 
riod 8: standard horseshoe (es(fc)) and horseshoe with a nega- 
tive full twist (es(b_i).) Orbits with the same exponent sum 
are braid conjugates. See [27] for the explicit conjugations. 
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TABLE III. Topological invariants for horseshoe braids up to period eleven (periods one through nine previously published 
by Hall [27]): name, code, permutation, Thurston type, rotation number, rotation interval, height, depth, exponent sum, and 
topological entropy of the braid calculated from a train track [11]. 
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